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Abstract 

We study the dynamics of a spin facilitated Ising model with long range kinetic constraints. To 
formulate those restrictions within an analytical approach we introduce the size of a kinetic active 
environment of a given spin. Based on a Master equation in second quantized form, the spin- 
autocorrelation function is calculated. It exhibits a pronounced slow dynamics, manifested by a 
logarithmic decay law of the spin-autocorrelation function. In case of an infinite kinetic interaction 
the mean field solution yields an asymptotic exact expression for the autocorrelation function 
which is in excellent agreement with Monte Carlo Simulations for finite interaction lengths. With 
increasing size of the active zone the cooperative processes, characterizing the facilitated model 
with short range kinetic interaction, become irrelevant. We demonstrate that the long range kinetic 
interaction dominates the actual spin configurations of the whole system and the mean field solution 
is the exact one. 
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I. INTRODUCTION 



A central topic in non equilibrium statistical physics is the quantitative and qualitative 
understanding of long-time phenomena in strongly interacting many-body systems without 
long range order. Such systems are characterized by a whole hierarchy of degrees of free- 
dom 1], which generate for example a broad distribution of relaxation times. The different 
levels relax in series resulting in an extreme slowing-down of dynamical processes. These 
processes and other singular behavior is in contrast to the critical slowing-down observed in 
the vicinity of conventional phase transitions. As a consequence of the mentioned hierarchy 
strongly disordered systems display a high degree of cooperativity of local processes. Thus, 
a signature for a broad class of dense systems are autocorrelation functions that can be 
fitted by a stretched exponential decay law for sufficiently long times or by a non-Arrhenius 
temperature dependence of relaxation times. Although it remains an open problem whether 
these phenomena are universal, they are, at least, a characteristic feature for many systems, 
in particular, when the systems are far from thermal equilibrium. For instance, such sce- 
narios arise for the main glass transition process in supercooled liquids. 
One way to model the hierarchy of relaxation processes is given by kinetic constraints. Due 
to that constraint a local dynamical process is free to evolve after changing its environ- 
ment. A rather simple realization is given by facilitated kinetic spin systems which offer the 
typical slowing-down of dynamical processes, caused by an increasing cooperativity of local 
spin-flip processes with decreasing temperature. In the present paper our interest is focused 
on the analysis the f-spin facilitated kinetic Ising models , originally introduced by 

Fredrickson and Andersen. These models are formulated on rf-dimensional lattices. To each 
lattice site i we assign a spin variable oi with two possible states Oi = ±1. The total set of 
all observables a = {o"j} forms a configuration. The basic dynamic of the facilitated kinetic 
Ising model is given by stochastic spin-flip processes <7j = +1 <-> <7j = —1. Additional to 
these stochastic Glauber dynamics p the system is subjected to self-induced topological 
restrictions. The origin of the imposed constraints becomes evident when one considers the 
local spin variable as coarse-grained local density variable. The down-state is associated 
to a low-density region while the up-state represents a high-density region. Whereas the lo- 
cal stochastic spin - flip processes are controlled by the thermodynamic Gibbs measure, the 
constraints reflect the mentioned cooperativity arising at low temperatures. An elementary 
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spin-flip process at lattice point i is only possible after an appropriate rearrangement of the 
spins in the neighborhood of the site i. In other words, a given spin <7j has to wait until the 
degrees of freedom in its environment have reached a certain configuration specified below. 
Insofar, the kinetic restrictions are related to the hierarchy of degrees of freedom (^J. 
To be specific, the restrictions are explicitly taken into account by choosing the local flip 
rates depending on the orientation of its neighbors. A spin at the lattice point % can only 
undergo a flip process in case the number z of nearest neighbored lattice point in an up 
state <Tj = +1 j_= 1 ... z is equal or exceed the restriction number / with < / < z. Such 
kind of models j2, S 0, Q are denoted as f-spin facilitated Ising models. For a realization 
on a d-dimensional lattice we use the abbreviation SFM[f, d]. Notice that in the original 
version the SFM[f, d] comprises a system of noninteracting Ising spins. A coupling is realized 
only via the kinetic constraints. Identifying the spin variables as local density variables as 
done above, the restrictions lead inevitably to cooperative rearrangements of the whole spin 
systems. Due to the kinetic constraints glassy effects arise at low temperatures (a, Isl LLO, LUj • 
The SFM[f,d] allows also another interpretation in terms of chemical reactions. To that aim 
let us identify the spin-up state as a molecule of type A and the down-state as a molecule of 
type B, respectively. Then the spin facilitated model can be considered as a reversible auto 
catalytic chemical reaction f A + B ^ (/ + 1)A on a d-dimensional lattice. The molecule 
changes its constitution from A to B and/or from B to A only in case, when at least / 
molecules of type A are situated in the nearest environment. 

Obviously the dynamics of the SFM[f, d] may not only determine by the restriction number 
/, but also by the range of the kinetic active environment of a certain spin <7;. In particular, 
the long time behavior of the SFM should be changed significantly, when the number of 
spins in the neighborhood is enlarged. One aim of the paper is to extend the original version 
of the spin facilitated Ising model in such a manner that a long range interaction will be 
included into the model. To realized a minimal extension of the SFM we define a quantity 
Ei(R) as the environment of an arbitrary spin Cj situated at the lattice point i. Here, the 
new parameter R is, so to say, the radius of the neighborhood. The simplest assumption 
consists of that Ei(R) contains all the spins cx^, the distance of which to the spin er^ fulfill 
the condition | k — % | < R. The total number of spins within that region of the radius R is 
denoted by G where G > f. Thus, if at least a subset of / spins in the region Ei(R) is in 
the up state +1, the spin <Tj is able to flip. As a consequence we expect a more pronounced 
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slowing down of the relaxation processes with an increasing ratio f/G with fixed / as well 
as for increasing / with fixed ratio f/G. Indeed, we demonstrate an utmost drastic slowing 
down, namely a logarithmic decay of the autocorrelation function. 



II. ANALYTICAL INVESTIGATIONS 
A. Evolution equations 

Let us start with a brief review of the main ideas behind the Fock-space formalism, 
which is a very powerful method for analyzing classical many-body systems with a stochas- 



tic dynamic given by a master equation on a lattice, for a review see [12j. The Fock-space 
approach is based on a quantum formalisms for the underlying master equation written in 
terms of creation and annihilation operators. This approach for SFM[/, d] will be more 
transparent by introducing occupation numbers via the relation <ji = 2rii — 1. Thus, the 
spin down orientation Oi = — 1 and up orientation oi — 1 correspond to an empty, rij = 
and a single occupied site with rii — 1, respectively. The SFM[/, d] can then be studied in 
the lattice gas representation with excluded volume effects, i.e. the change of a certain con- 
: iguration n = in,- 1 is only possible when the exclusion principle is fulfilled. Following Refs. 
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I15I 16 , [17I, [l^l] , the probability distribution P(n, t) can be related to a state vector 
\F{t)) in a Fock-space according to P(n,t) = (n\F(t)) and \F(t)) = J2 n P(. n ^) l n )' respec- 
tively, where the basic vectors |n) can be expressed by second quantized operators. Using 
this representation, the underlying master equation can be transformed into an equivalent 
evolution equation in a Fock-space, written in the form: 

d t \F(t)) = L\F(t)) . (1) 

The dynamical matrix L(n, n') of the master equation is mapped onto the operator L = 
L(a, a^), where a and a* are the annihilation and creation operators, respectively. Originally, 
his transformation had been applied for the Bosc case with unrestricted occupation numbers 

inn nnn 

1J,114|,113- Here, we consider the case of restricted occupation numbers |16l. 1171. 1181] . In order 
to preserve the restriction of the occupation number in the underlying dynamical equations 



the commutation rules of the operators a and a) are chosen as Pauli-operators 
with the commutation rules: 
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[di, a]] = Sij(l - 2a\ai) [a { , aj] = [a], a\] =0 a? = (aj) 2 = . (2) 
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The master equation of the SFM[/, d] can be expressed by the following evolution operator 
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) ^\k 1 ...k f N kl N k2 ...N kf 1 \p{oi - Ni) + A(a} - (1 - JV<))] , (3) 



'GV 1 

ki,...kf 

where iV; = ajaj is the particle number operator which satisfies iVj|n) = nj|n). Whereas the 
temperature dependent flip-rate A stands for the process |— >t, the rate (3 > A moderates 
the inverse process. The temperature T is introduced in the model by T ~ [ln(/3/A)] _1 . 
The terms in the squared bracket in Eq. (jSJ) represent a single spin-flip process at lattice 
site i. The product of the particle number operators N kl Nk 2 ...Nk f takes into account the 
constraints for the flip-process at site i. The lattice function Ki\ kl ...k f is only different from 
zero if all / lattice sites ki . . . kf are different from each other and moreover, all sites are 
inside the kinetic active zone around the lattice point i, i.e. all sites k a with a = 1 . . . / 
satisfy the condition \k a — i\ < R. In that manner flip processes are only allowed if all 
the lattice sites k a , a — 1 . . . / in the surroundings of site i are occupied. Notice that the 
restriction number / with / < G, where G is the number of considered neighborhood sites 
in the environment of the lattice point i. Introducing the operator 

S= Yl N k 

keEi(R) 

the constraints can be rewritten using some combinatorial manipulations. Then the evolu- 
tion operator Eq. Q reads 

L = E (/) 1 (/) - N <) + A ^ f - o- - N ^} ■ ( 4 ) 

Obviously, the rate for the flip-process a = +1 ^ a = — 1 becomes (3 for the case that all 
spins are in the up state, independently from the choice of Ei(R) and therefore of G. To 
proceed further we follow Doi [l^J and calculate the average of an arbitrary physical quantity 
B(n) by using the average of the corresponding operator B{t) = J2 n l n ) -^( n ) ( n l via 1^3 1 

(B(t)) = ^P{n,t)B(n) = (s\B\F{t)). (5) 

n 

Here we have used the reference state (s\ = (nl. The normalization condition for the 

n 

probability density is included in the condition (s\F(t)) = 1 with the consequence [23], 



that the evolution operator fulfills always the relation (s \ L = 0. Note, the special feature 
of this Fock-space formulation is due to the fact that the average value is linear in the 
corresponding state vector in contrast to quantum mechanics where it is bilinear. In the 
same way, correlation functions can be expressed as 

(A(t)B(t')} = ^A(n)P(n,t;n',t')£(n') = (s \ Aexp {L (t — t')} B\ F(t')) . (6) 

n,n' 

Furthermore, because of Eq. and Eq. (0) the evolution equation for an arbitrary operator 
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B(t) } for example the particle number operator, is given by| 

d t (B) = (s\[B,L]\F(t)) , (7) 

which can be extended immediately in order to get the kinetic equations for time-dependent 
correlation functions. As a general result of the quantum formulation of stochastic processes, 
all dynamical equations, describing the classical problem, are completely determined by the 
commutation rules of the underlying operators and the structure of the evolution operator 
L. The Fock-space approach for the master equation allows t 

nr 

of evolution processes, e.g. aggregation, chemical reactions |17Lllj 

as well as the spin facilitated kinetic Ising models, compare also [12| . One of the advantages 
of the Fock-space approach is given by the simple construction pri nciples for each evolution 
operator L on the basis of creation and annihilation operators |25j |. 
The evolution equation for the averaged particle number operator, using Eq. (|7J) reads 



re analysis of a broad variety 
non-linear diffusion 124 



-i 



s\\ /fs\,\\ / fs 



(8) 



The evolution equation for the autocorrelation function becomes 

(A Vf^WwwV ( 9 ) 



.// w / 

which is valid for t > t' . Due to the constraints both evolution equations contain higher 
moments of spins at neighboring lattice sites. The evolution equation of these moments 
contains the next higher correlations, i.e. one obtains an infinitely large hierarchy of evolu- 
tion equations. In Eq. (jHJ) the higher moments decouple in the equilibrium state, because 
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the Hamiltonian of the SFM[/, d] is that of a simple paramagnetic gas. Furthermore, the 
equilibrium is given by d (N{) jdt — leading immediately to 



This simple decoupling fails outside the equilibrium due to the strong kinetic coupling be- 
tween neighbored spins. But in case of diverging size R of the active zone and consequently 
for G — > oo the intrinsic self averaging effect due to the infinite number of neighbored spins 
results again in a decoupling of higher correlations. 



B. Long range kinetic interaction 

When the system is subjected to homogeneous initial conditions then the averaged particle 
number operator is independent on the lattice site. Denoting that value by (iVj) = iV we 
obtain for large G (G — > oo) from Eq. (JHJ) 



" ( G / iV ) (A [1 - N] - 0N) — (A + P) g) " ( G ^) (N„ -TV) . (10) 
Using the same arguments the evolution equation for autocorrelation function, Eq. (JHJ), reads 



ON _ / G 



| (N(t)N(t')) = f f) 1 ( GN f (t) ) (A [iV(t') - (N(t)N(t'))} - P (N(t)N(t'))) . (11) 



The last equation can be simplified by introducing the two-time correlation function 



H{t,tf) = N(t)N(t') - N(t)N(t'). It follows 



d .„ , „ x fG\' L [GN(t)\ 



di H(t 1 lf) = -(X + P)l \ ( f KJ )H(t,t'). (12) 



This equation can be solved explicitly. We get for t > t' 

, m m /^V 1 /" ^iV(r) . , 



Using Eq. (JTUJ), we find 



JV(t) _ _ 

1° 777T~~T = — (A + 3) / ^ _ = in A - (A + « 



_7 A-(A + /3)iV A - (A + /?) JV(f ) 

JV(t') 
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which yields 



H{t,t') N(t)-N ( 



cq 



H(t',t') N(t' 



N, 



(13) 



cq 



The correlation function H(t, t') depends on t and t' with t' < t. Whereas t is the observation 
time, t' can be interpreted as waiting time tw- The system relaxes from the initial state at 
t = into the instantaneous state at t via an intermediate time t' = tw- It seems to be 
reasonable to introduce the reduced correlation function C(r,tw) by the definition 

H(t w + r, t w ) 



C(r,t 



w , 



(14) 



H(t w , tw 

The reduced correlation function describes the relaxation behavior of the SFM[/, d]. 
The time dependent decay N(t) follows from the solution of Eq. (JTHJ). We start from the 
initial condition <Ji = +1 for all lattice sites i, i.e. we have the initial condition N(0) = 1 > 
A^q. Furthermore, one concludes from Eq. (jl()() the validity of 1 > N(t) > N eq for all times. 
The SFM[/, d] approaches the thermodynamic equilibrium state, if GN eq > f, otherwise the 
system reaches the stationary state N st = f/G > N eq . 



III. RESULTS 



A. Relaxation into the equilibrium state 



Let us firstly analyze the behavior for finite /. In this case we obtain for f/G — ► 0, i.e. 
G — ► oo and consequently 

] ' ] ^N. (15) 



./ 



/ 



Inserting this expression in Eq. (jlOj) the evolution equation reads now (compare also 2l| ) 



dt 

From here, we obtain the implicit solution 

'N-N ( 



N 



N, 



cq 



1 - 



A' 



N, 



C<| 



(16) 



-In 



' cq 



l-N, 



cq 



/ ]\[k-i r 



k=2 



1 - 



N 



k-l 



Ni q (\+p)t. 



(17) 



The sum can be rewritten by using the Lerch-$-function 

®(z,a,b) = 



n=Q 



(b + n) a ' 



(18) 
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in the form as follows 

^$ (^p, 1, fj = $(iV eq , lJ) + (X + 0)t. (19) 

The function N(t) decreases monotonously from N(0) = 1 to N(oo) = N eq . The short time 
behavior can be calculated by a series expansion of Eq. (fTTj) . Up to second order in t we get 
for t -> 

N=l-pt+(f/3 + \ + (3)(3t 2 /2 + o(t 3 ) . 

This short time regime is valid for t < The early stage, described by the first 

order of the expansion, depends only on the kinetic constant quantifying downward spin 
flips. Backward flips and kinetic restrictions becomes relevant at later stages, because the 
corresponding parameter / and the kinetic rate appear firstly in the second order term. 
The final time regime is given by N N eq . Obviously, the first logarithmic term in Eq. (fTTjl 
is the dominate one. Then we find the asymptotic solution 

N(t) - N cq = (1 - N cq ) exp {- (X + 0) Ni q t} . (20) 

From Eq. (j!7|) follows, that this exponential decay regime becomes relevant for t ^> 

Further, the evolution equation allows a pronounced intermediate regime resulting for 
N eq <C 1 and large /. In this case we find for 1 ^> N(t) ^> N cq a power law behavior 
in form of 

li + fpt]~ 1/f ~r 1/f . (21) 

For infinite large /, but f/G — > both, the description of the decay of the initial regime 
and the intermediate regime may be combined to the general relation 

N=l-j\n[l + f(3t] (22) 

for all finite times This result can be obtained directly from Eq. ()16|) considering 

N(t) ss 1. 

In case of a finite ratio f/G and G — > oo we get for N f /G instead of (j!5)l the relation 

which indicates a further slowing-down of the relaxation with decreasing ratio f/G and 
increasing N. An explicit expression for the time evolution can be derived by inserting 
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Eq. J2ED in Eq. (fl) and making the ansatz N(t) = 1 - y(t) with y < /3/(A + /3) < 1. It 
results 

2G \Pf 
f 2 " 



iV(t) ~ 1 - — In 



-t 



(24) 

e R reauires nu 



2G 

While the case of R —>■ oo can be solved analytically, the analysis for finite R requires numer- 
ical methods. Here, we use a standard Monte-Carlo Simulation 
dimensional lattice. Fig. ^ shows the decay N(t) for fixed G and various / <C G in compar- 
ison with the predictions obtained from the mean field theory. As expected, the accuracy of 
the mean-field solution Eq. increases with decreasing ratio f/G. This phenomenon is 
partially justified by the relation Eq. (|2*H|) which slows down the decay described by Eq. (|19j). 
Additionally, an increasing ratio f/G requires an increasing number of spin flips to collect / 
up-spins in the environment of the lattice site % before the spin Oi is able to flip. While the 
first effect is contained in the general mean field approach given by Eq. (|10p. the cooperative 
phenomenon is not taken into account by the mean field concept. Fig. El shows the decay 
of N(t) for various / and G = If . Because the temperature T ~ (ln/3/A) -1 is very small, 
the initial regime is partially comparable with the situation for T = 0. Here the dynamics 
stop at a non ergodic state. The approach into the non ergodic state can be demonstrated 
by rigorous results j^J, in this connection see also 2(|. In particular, for / = 1 and d = 1 
the decay of the ordered phase N = 1 is given by 

N (t) = exp {exp {-(3t) - 1} . (25) 

In the long time limit this results in a pronounced slowing-down leading to a non ergodic 
behavior which is manifested in that for t — ► oo the quantity N(t) remains finite with 
N(oo) = e' 1 . 

For all finite temperatures the non ergodic state is unstable and decays by cooperative re- 
arrangements of the spin configurations. For sufficiently long times the system approaches 
the equilibrium state, N(t — > oo) = N eq . This slow process is well described by a mapping 
of the kinetically constrained spin flips onto an auto catalytic reaction diffusion process on 
a one-dimensional lattice (3], which offers a slightly stretched exponential behavior. 
The decay N(t) approaches again the mean field result, obtained in Eq. (flH|l for sufficiently 
large /, which indicates the decreasing relevance of cooperative processes for large /. How- 
ever, a larger difference between the solution Eq. ([19)1 and the numerical result remains. 
This difference occurs for all N(t) < 1 and it is related to the fact that Eq. (|T5jl is valid 
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only for f/G — > 0. For finite f/G, Eq. (JT5|l is an approximation which becomes exact only 
for N(t) — > 1, see Eq. (j^SJ. Therefore, the decay limit Eq. (J2"2"|) obtained for / — » oo is also 
valid for all finite values f/G with t <C e^. 



B. Correlation function 

The knowledge of N(t) allows us to compute the reduced correlation function C(r,tw), 
defined by Eq. ()14j) . The decay of this correlation function is also characterized by different 
regimes. 

(i) For a short waiting time tw — > and short time scale r we find an universal relation for 
the reduced correlation function, namely 

C(r, t w ) = 1 - (A + p) (1 - (5ft w )r . (26) 

Note, that a general analysis of Eq. (JSJ) and Eq. (JSJ) shows that this quantity is independent 
of the interaction range R. 

(ii) The limit of a long waiting time tw — > oo leads always to 

C(r,t w )=exp{-(\ + (3)Ni q r} . (27) 

The independence of tw indicates the proximity to the thermodynamic equilibrium. 

(iii) If the time scale increases, i.e. r — > oo, we obtain 

C(r, t w ) = exp {- (A + /5) iV/ q (r + t* w )} (28) 

with 

ln(iV(%) - iV eq ) 

(A + ^)iV/ q • 

(iv) Within a short or intermediate time scale r and an intermediate waiting time t w the 
reduced correlation function offers a characteristic aging behavior manifested by 



C(r,t 



Wi 



t w 



-1// 



(29) 



The crossover to infinite / shows, that the decay processes freeze in. This means, the short 
time regime is restricted in terms of tw, t < (/9/) _1 . In the intermediate time regime the 
correlation function becomes 
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The result is valid up to r/tw ~ e ■ The correlation function C(r,tw), characterizing the 
aging behavior, reveals a logarithmic decay with an extreme small prefactor. Thus, the 
intermediate time regime dominates the behavior of the system for all finite times r and tw- 
Our numerical investigations support these mean field results. Fig. |3] shows various corre- 
lation functions directly obtained from the numerical simulations in comparison with the 
correlation functions Eq. (jl4j) computed from the numerically determined decays N(t) using 
Eq. (j!3|) . As discussed above in the context of the numerical analysis of the decay function, 
significant differences appear only for small /. 

Fig.0]shows the dependence of the correlation function on the waiting time tw for f/G = 1/2 
and / = 50. We obtain a typical crossover from an exponential decay, independent on the 
waiting time to a power law like behavior which scales with the waiting times tw- This 
behavior agrees again with the mean field approach. 

Finally, we demonstrate the dependence of the correlation function on the restriction number 
/ for a fixed waiting time, see Fig. 03 As discussed above, we observe a continuous crossover 
to a logarithmic decay for sufficiently large /, suggested by Eq. (fHUJ) . 

IV. CONCLUSIONS 

We have studied the influence of the restriction number / and the range of the nearest 
environment on the kinetics of the SFM[/, d]. To this aim the original model was extended 
by an additional parameter G describing the number of kinetically neighbored spins around 
a given lattice site i. All spins of this environment have the same influence on the flip prob- 
ability of the spin <jj. We obtain a pronounced slowing down of the relaxation behavior and 
the correlation functions for an increasing ratio f/G and fixed / as well as for increasing / 
and fixed ratio f/G. For sufficiently large /, a logarithmic decay of the spin-spin correla- 
tion function was observed. Furthermore, in the case of an infinite kinetic interaction range 
G, the correlation functions in the kinetic equations decouple and the mean field approach 
becomes an exact theory. This approach allows the identification of various regimes. While, 
the initial regime is described by universal relations which are valid also if the decoupling 
procedure fails, the final regime is always an exponential decay describing the relaxation 
into the equilibrium state. The intermediate regime can be approached by power laws for 
the correlation functions, which are significantly controlled by the waiting times. 
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For finite G we observe a further slowing down due to the fact, that cooperative processes 
become relevant. These cooperative rearrangements are necessary to collect / up-spins in 
the environment of the lattice site i so that the spin o"j is able to flip. This effect appears 
especially for small kinetic ranges G and large ratios f/G. In the case of increasing G the 
kinetically active environment E± of a given lattice site % becomes more and more represen- 
tative for the actual spin configuration of the whole system and the mean field solutions 
converge to exact results. 
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FIG. 1: N(t) for various / and fixed G = 200 in comparison with the analytical mean field results 
(straight lines). The temperature is given by T = (3 In 10) _1 . The values of / are / = I (squares), 
/ = 2 (up triangles), / = 5 (circles), / = 10 (down triangles), / = 20 (crosses) and / = 50 
(diamonds) 
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FIG. 2: N(t) for various / and fixed ratio f/G = 1/2 in comparison with the analytical mean field 
results (straight lines). The temperature is given by T = (31nl0) _1 . The values of / are / = 1 
(squares), / = 2 (up triangles), / = 5 (circles), / = 10 (down triangles), / = 20 (diamonds) and 
/ = 50 (crosses) 
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FIG. 3: Correlation function C(tjy> r ) computed from Eq. (j!3[) using the numerically determined 
function N (straight lines) and obtained from direct numerical simulations for tw = 1 an d t\y = 10 6 
(/ = 10, G = 200: squares and circles), (/ = 10, G = 20: up and down triangles), (/ = 1, G = 2: 
crosses and diamonds) 
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FIG. 4: Correlation function C(t\y,T~) for fixed / = 50 and G = 100 and different waiting times 
t w = 10™ with n = . . . 6. 



17 




18 



